
  #include "stdio.h"
  #include "5strq.c"
  main()
  { int i,j;
    double b[5],c[5],q[5][5];
    double a[5][5]={ {10.0,1.0,2.0,3.0,4.0},
        {1.0,9.0,-1.0,2.0,-3.0},{2.0,-1.0,7.0,3.0,-5.0},
        {3.0,2.0,3.0,12.0,-1.0},{4.0,-3.0,-5.0,-1.0,15.0}};
    strq(a,5,q,b,c);
    printf("MAT A IS:\n");
    for (i=0; i<=4; i++)
      { for (j=0; j<=4; j++)
          printf("%13.6e ",a[i][j]);
        printf("\n");
      }
    printf("\n");
    printf("MAT Q IS:\n");
    for (i=0; i<=4; i++)
      { for (j=0; j<=4; j++)
          printf("%13.6e ",q[i][j]);
        printf("\n");
      }
    printf("\n");
    printf("MAT B IS:\n");
    for (i=0; i<=4; i++)
      printf("%13.6e ",b[i]);
    printf("\n\n");
    printf("MAT C IS:\n");
    for (i=0; i<=4; i++)
      printf("%13.6e ",c[i]);
    printf("\n\n");
  }

